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Abstract 

We investigate the magnetic moment operator for constituent fermion masses 
for chirally symmetric theories. Constituent fermion masses are generated through 
a yukawa interaction of the fermion with a scalar (or /and psuedoscalar) field via 
the vacuum expectation value (VEV) of the scalar (or and psuedoscalar) field. 

We especially consider the high baryon density ttq condensed phase, in which 
chiral symmetry is spontaneously broken, with space varying expectation values 
of the a and ttq fields. This phase has a spin polarized fermi sea as the ground 
state. We show that there is indeed generated a macroscopic magnetization in 
this phase, contrary to what one would have found, if one just used a primitive 
phenomenological magnetic moment formula for explicit/ current fermion masses. 

Furthermore, this analysis reveals that the magnetization of this state goes up as 
the VEV, that determines the 'mass', comes down with increasing baryon density. 
The consequent high magnetic field that is generated will destabilize this state at a 
threshold density. This is important in the context of neutron stars, as such a high 
density state may be responsible for very high magnetic fields in the dense core of 
neutron stars. This could potentially be the origin of magnetars - the stars with the 
largest magnetic fields in the universe. 
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1 Introduction 



Constituent masses are generated via yukawa interactions of the fermion with a scalar (or 
/and psuedoscalar) field via the vacuum expectation value (VEV) of the scalar (or and 
psuedoscalar) field. Constituent fermions arise in a variety of models from the Standard 
model of the electroweak interactions via the Higgs VEV or in the strong interactions(in 
the chiral symmetry limit) where the nucleon / quark masses are generated by the VEV 
of the sigma- field < a >. Constituent masses are different from explicit or current masses. 
The obvious difference is that under a chiral symmetry transformation the exact chiral 
symmetric Hamiltonian remains invariant (constituent masses) whereas any current mass 
in the Hamiltonian explicitly breaks the chiral symmetry. 

Also, there are physical situations in which the VEV's of the scalar (or and psue- 
doscalar) field and therefore constituent masses can vary, like at finite density or tempera- 
ture, which underlines the difference with current masses which 'do not' change with finite 
density or temperature. Indeed, there are situations where the VEV's can even become 
space varying, making the difference with current masses even more dramatic. In such a 
context the use of the usual phenomenological formulae(e.g. for magnetic moment), that 
do not take account of this, can be misleading. 

One such case that we will consider is the case of a high baryon density ground state 
for strongly interacting nucleons, which has a ttq condensate that is a (space varying) sta- 
tionary wave | Dautry and Nyman 1979|| Baym 1977|Kutschera, Broniowski and Kotlorz 1990 



Soni and Bhattacharya 2006] . In this case when we calculate the magnetic moment of 



the nucleons in the presence of a ttq condensate using the naive formula, we find a some- 
what strange result: that for a spinpolarised neutron ground state, the magnetic moment 
vanishes when averaged over spacd_|. However, when we compute the magnetic moment 
operator from first principles, from the chiral lagrangian, we find it has the right space 
dependence, which cancels out with the space dependence of the ground state to give a 
magnetic moment that is proportional to the total spin. Furthermore, this reveals that 
the magnetization of this state goes up as the VEV, that determines the 'mass', comes 
down with increasing baryon density. 

This is important in the context of neutron stars, as such a high density ground 
state may be responsible for very high magnetic fields in the dense core and could be the 
origin of magnetars | Bhattacharya and Soni 2007lH"aridass fc Soni 20TU] - the stars with 



the largest magnetic fields in the universe. 
1 ~R. F. Sawyer private communication 
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2 Current and Constituent Masses 



2.1 QED 

Quantum Electrodynamics is a theory with an explicit or current electron mass, m e , which 
breaks chiral symmetry explicitly. The case of the electron magnetic moment operator 
is explicitly worked out in the text of Sakurai (Advanced Quantum Mechanics) via the 
Gordon decomposition (see (3-3 page 85) | Sakurai] ). 



el 

H mag = ~^^F ^ . (1) 



(h 



2.2 The Chiral symmetric Gellmann Levy sigma model 

In a chiral symmetric theory, the mass of the nucleon/quark comes from the VEV of 
< a >, for example in the linear cr-model of Gellman and Levy[Gellmann Levy I960]. In 
the case of exact chiral symmetry we have the following lagrangian when we couple the 
Gellmann Levy sigma model to the electromagnetic field 



L = $ {Jp + g y (a + i l5 r • tt)) ^-\{d^a) 2 -\{d^) 2 -^{a 2 +n 2 -{F) 2 ) 2 

(2) 

In the limiting case of a small explicit pion mass being set to zero, the usual (uniform 
in space) symmetry breaking that follows on the minimization of the potentials above is, 
< a >= F = /^and < n >= 0. 

The masses of the scalar (PS) and fermions are given by 

< o > 2 = F 2 (3) 
where F is the pion decay constant. It follows that 

ml = 2X 2 F 2 m = g < a >= gF (4) 

The quark or nucleon mass m, is a spontaneous mass that is generated from the 
spatially uniform VEV. In this case the usual magnetic moment formula above ([T]) works. 

On the other hand when the VEV's of < a > and < tx > depend on space coordi- 
nates, the above expression for the magnetic moment will not work as we will find in the 
following section. 
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3 The 7Tq condensation : space dependent VEV's 



3.1 The ttq condensed ground state 

Here we shall consider another realization of the expectation value of < o > and < tx > 
corresponding to ttq condensation. This phenomenon was first considered in the context 
of nuclear matter | Dautry and Nyman 1979 , Baym 1977| . Such a phenomenon also oc- 



curs with our quark based chiral a model and was considered at the mean field level by 
Kutschera, Broniowski and Kotlorz for the 2 flavour case |Kutschera, Broniowski and Kotlorz 1990| 



and by one of us for the 3 flavour case Soni and Bhattacharya 2006 . Working in the 



chiral limit they found that the pion condensed state has lower energy than the uniform 
symmetry breaking state we have just considered for all density. This is expected, as the 
ansatz for the pion condenced phase is more general. 

The expectation values now carry a particular space dependence 

< a > = F cos (q.r) (5) 

< 7r 3 > = F sin (q.r) (6) 

< 7Ti > = (7) 
<vr 2 > = (8) 

Note that the relation, < a 2 > + < tt 2 >= F 2 , is preserved under this pattern of 
symmetry breaking. Also, when \q\ goes to zero, we recover the usual space uniform phase 
above. 

The Dirac Equation in this background is solved in | Dautry and Nyman 1979 



Kutschera, Broniowski and Kotlorz 1990| by the artifact of writing the wavefunction, ip, 



in terms of a chirally rotated wave function, x(k), 

i/j = exp(-i(r 3 /2)7 5 g • r) ■ X (k) (9) 

where, X (k) are momentum eigenfunctions. 
The Hamiltonian reduces to 

H X {k) = (a.k - ^q.a l5 T 3 + 0m) X {k) = E(k) X {k) (10) 



where m = g\/ < a > 2 + < 7r > 2 = gF 

The second term arises from the condensate and has been written in terms of the 
relativistic spin operator, cr/5 . It is evident that if spin is parallel to q and r 3 = +1 
(proton/up quark) this term is negative and if r 3 = —1 (neutron/down quark) it is 
positive. For spin antiparallel to q the signs of this term for r 3 = +1 and —1 are reversed. 

The spectrum for the hamiltonian is the quasi particle spectrum and can be found 
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to be | Dautry and Nyman 1979 



Kutschera, Broniowski and Kotlorz 1990 



E H {k) = Jm 2 + k 2 + ^q 2 - ^ m 2 q 2 + (q.k) 2 (11) 

E {+) (k) = ^m 2 + k 2 + ^q 2 + \Jm 2 q 2 + {q.k) 2 (12) 

The lower energy eigenvalue Er^ has spin along q for r 3 = 1, or has spin opposite to q 
for r 3 = — 1 . The higher energy eigenvalue has spin along q for r 3 = — 1 , or has spin 
opposite to q for r 3 = +1. 

The ground state is constructed by occupying all the lower energy states Et \ as the 
gap between the lower and higher states is large. In this background the fermi sea is no 
longer degenerate in spin but gets polarized into the states above. 

3.2 The naive magnetic moment calculation in the ttq condensed 
state 

Let the condensate wave vector q, define the z axis. For neutrons r 3 = —1. Since the 
spins are aligned opposite to q, we expect that there should be a net magnetic moment 
opposite to the q axis . Let us look at the source of the magnetic field, from the 'usual' 
Gordon decomposition of the electro magnetic current j Sakurai] . 

g 

H m ag = ^~ F ^ ' ^ (^,uT 3 ) ^ (13) 



where, m = g\J < a > 2 + < 7r > 2 . 

Substituting for ip = exp(— 2^3/2)75^ • r) ■ x(k), where, x(k) are momentum eigen- 
f unctions, we find 

g 

Hmag = 7^ F >^ ' x(k)a^ u T 3 exp(-zr 3 7 5 g ■ f)x(k) (14) 

The ground state is constructed by filling the spin polarized E_(k) states up to the 
fermi momentum. The total spin of the ground state is the sum over the spin polarized 
Fermi-sea. Do the magnetic moments also add up like the spins? 

If we simply used the normal (explicit current mass) formula fTT4|) we would get a 
space average over the above sing ■ x and cosg ■ x terms that goes to zero. 

3.3 The full chiral magnetic moment calculation in the ttq con- 
densed state 

Using the notation of Sakurai | Sakurai] . 7 M <9 M = 7i<9j + 74^4 = 7^ + (ijo)(— id ) we may 
write the following chiral Dirac equations for the the ttq condensed state 
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(l^dy - ie^/yAy + g(a + ij 5 r -7f))V> = (15) 
(7m^ _ iej^An + mexp(ir 3 7 5 g • r))^p = (16) 

Now we can invert this equation to get 

■ip = — exp(-zr 3 7 5 g ■ r)^^ - iej^A^i/; (17) 
Similarly, by complex conjugation with a little algebra we find 

rp = —{dpi/ryu + ieifrypAp)) exp(-ir 3 -f 5 q ■ f) (18) 

Following the 'Sakurai' trick j Sakurai] . we use the ploy of dividing the electro mag- 
netic current operator into two identical parts 

g 

Jm=2^ 7 ^ + ^ 7 ^ ( 19 ) 

We then use the chiral Dirac equation fll7p above for if), and insert it in the first term 
and (I18p for , ip, and insert it in the second term. 

After doing the Gordon decomposition, we get the following expression for the elec- 
tromagnetic current operator (if we leave out all the terms that depend on the vector 
potential, Ay), 

g 

3n = 2~ exp(2r 3 7 5 g • r)a^ u T 3 ip + ip exp(ir 3 7 5 g • r)a^ u T 3 (d u ip)} 

+ exp(zr 3 7 5 g • r)r 3 ip - ip exp(zr 3 7 5 g • r)T 3 (d^))) (20) 

Notice the presence of the phase factor in the current operator. The expression 
reduces to the usual operator, if we put q = ,for then the phase factor drops out (= 1). 
The first part has the Lorentz structure of the magnetic moment operator and the second 
term reduces to a current density in the Schroedinger theory [ Sakurai] . For our present 
purpose we will concentrate on the first term which is the magnetic moment type term. 

g 

Hmag = -j™ 9 Ay = ~i^Ay [d u (i/j exp {it 3 7 5 <f ■ r)ay iU T 3 ip) 

- ^(igi/TsTa) exp(zr 3 7 5 g ■ r)a y, tV T 3 ijj\ (21) 

By partial integration of the first term above, we can move the derivative to act 
on the vector potential , Ay, and write it in the usual form of the magnetic moment 
term. However, now there is an extra term, the second term, which needs a physical 
interpretation. 
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The first term 

The first term is the regular magnetic moment term and is 

g 

Himag = —F^ v [ip exp(iT 3 j 5 q ■ r)a^ u r 3 ij] (22) 

On substituting for if) = exp(—i(r 3 /2)^f 5 q-r)-x(k) we find that the two phase factors, 
from the magnetic moment operator and the wavefunction, nicely cancel to give 



— ■ X (k)a^x(k) (23) 

This magnetic moment is not space dependent and does not average to zero and 
thus restores the proportionality between the total spin and total magnetic moment for 
the ground state. 

The second term 



H 2 mag = -(—)A fi [-i/;(iq u j 5 T 3 ) exp(ir 3 7 5 g • r)cr Mii ,r 3 ^] (24) 

By using the field equations above for ip and ip and the 'Sakurai trick' | Sakurai"] the 
second term can be recast to yield many terms. 

We drop the terms that are not relevant for an effective magnetic moment term for 
reasons of economy. We also drop terms which are quadratic or higher order in 'q' and / or 
A^. We are then left with only a single term, which has the form of a magnetic moment 
operator. Below we write this term in terms of X'- 



(i %^) F ^ ' x(k)q u i 5 T 3 a^ x x(k) (25) 

This term is a somewhat different Dirac bilinear that interacts with the magnetic 
field. We note in passing that the electromagnetic interaction is not invariant under chiral 
transformations. 



4 The magnetic field in the tto condensation phase 
transition as a function of baryon density 

The interesting question that this poses in the context of the ttq condensation phase 
transition is how does the chiral mass term behave as a function of baryon density? It 



can be seen from fig 4 Kutschera, Broniowski and Kotlorz 1990 that as the density goes 
up the value of the chiral mass, m = g\/ < a > 2 + < 7? > 2 , decreases monotonically. 
As this mass appears in the denominator of the expression for the magnetic moment, 
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the magnetic moment keeps going up as the density increases. This implies that the ir 
condensed state in which spins and magnetic moments are aligned would have higher 
magnetic moments and therefore higher magnetic fields as the density goes up. 

Will the effect keep increasing as density goes up or is there a cap on the magnetic 
field? As baryon density is increased the chiral mass parameter goes down, lowering the 
condensate energy but simultaneously increasing the magnetic energy. An upper bound to 
the magnitude of the allowed magnetic field can be obtained by equating the condensation 
energy density to the magnetic energy density 

777 

(A 2 /4) • (-) 4 * B 2 /(8n) (26) 

Kutschera et al |Kutschera, Broniowski and Kotlorz 1990| use a sigma particle mass 
of 1200M ev corresponding to A 2 ~ 75. For a more realistic sigma particle mass of QOOMev 
we have A 2 ~ 20. From the figure 4 in |Kutschera, Broniowski and Kotlorz 1990 the 



value of the chiral mass parameter at a baryon density of ~ 1/F 3 is about half its value 
at density (m = gf^). This yields an upper bound for the magnetic field that the 
condensate can sustain, ~ 10 17 gauss, which is close to the field that we found ~ 10 16 gauss 
| Bhattacharya and Soni 2007 lHaridass fc Soni 20T0]. as the core field in our magnetar 



model. The implication is that in the tiq condensed state in the core of the star we are 
close to upper critical field. Beyond this density /magnetic field the ground state is not 
stable and is likely to transform into a charged pion condensate or unpolarized fermi sea 
( which carry no magnetic moment) - an interesting inference. 
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